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Abstract 

In the present study, we use a fifth-order ordinary differential equation, a generalization of the Blasius 
equation derived from the Leray-a model, to examine the transition to turbulence in rough-wall boundary- 
layer flows. This equation, combined with a weaker formulation of the von Karman log law modified 
to include the effects of surface roughness, provides a family of turbulent velocity profiles with two free 
parameters: the momentum thickness based Reynolds number and the roughness function. As a result, 
a skin-friction correlation is obtained and predictions of the transitional Reynolds numbers and maximal 
skin-friction values are made based on the roughness function. 

PACS numbers: 47.27. -i, 47.27.Cn, 47.27.E-, 47.27.em, 47.27.nb 
Keywords: turbulent boundary layer, rough wall, Leray-alpha model 



*Electronic address: lache skid @ umich . edul 
^Electronic address: ldkma@urmch.edul 



1 



I. INTRODUCTION 



Wall-bounded turbulent flows continue to be significant in both the natural environment and 
in engineering applications. As such, the turbulence community has persistently endeavored to 
describe the fluid dynamics in such flows ll^ . In particular, the determination of skin-friction (or 
wall shear stress) has been of interest to many researchers. 

For smooth- wall turbulent boundary layers, both direct and indirect techniques have been de- 
veloped to determine skin-friction. These include the use of the momentum integral equation, 
correlations based on pressure measurements at the surface, and fitting mean velocity profiles 
based on a defect or power law 



However, in the rough- wall case [fioi lllll. direct measurement of skin-friction is often difficult. 
Indirect methods have included using friction velocity, Ur, to estimate skin-friction. The modified 
Clauser method, for instance, has been used to approximate by fitting a logarithmic velocity 
profile to experimental data. However, as discussed by Acharya and Escudier yj], this technique 
is subject to large uncertainties because the degrees of freedom are increased from one (c/) for a 
smooth surface to three (c/, e, Au/ur), where c/ is the skin-friction coefficient, e is the error in 
origin, and Au/ur is the roughness function. Alternative indirect techniques include determining 
Ur using a velocity defect law, or power law formulations. Other correlations such as Bergstrom 
et al.'s skin-friction correlation with the ratio of the displacement and boundary-layer thicknesses 
have also been suggested [|2D. 

In llvl], a new theoretical method to derive a skin-friction correlation for smooth- wall turbulent 
boundary layers was presented. In this note we adopt it to rough- wall boundary-layer flows using 
the boundary-layer approximation of the Leray-a model of turbulence i4\ and the von Karman 



log law for rough walls [j8|, IE [l2|]. The benefit of this approach is that it leads to a prediction of 
the critical Reynolds number based on momentum thickness -R^"*, the minimal Reynolds number 
where the transition to turbulence may occur. More precisely, we obtain the following dependence: 

Rf^ = -51.8Au/ur + 365.5. 

Since in a turbulent boundary layer, the skin-friction coefficients attains its maximum cj^^ at 
Re = Rq^^^, we also obtain the the following relation: 
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II. BOUNDARY-LAYER APPROXIMATION OF THE LERAY-a MODEL 



Proposed as a closure scheme for the Reynolds equations [|4|], the Leray-a model is written as 

d 

—V + (u ■ V)v = i^Av -Vp + f 

V ■ u = (1) 

V = u — a^Au, 

where u is the averaged physical velocity of the flow, p is the averaged pressure, / is a force, and 
> is the viscosity. The filter length scale a represents the averaged size of the Lagrangian fluc- 
tuations and is considered as a parameter of the flow. More specifically, we assume that a changes 
along the streamlines in the boundary layer, and is proportional to the thickness of the boundary 
layer (see [|50). Inspired by the Navier-Stokes-a model (see yO and references therein), this model 
compared successfully with experimental data from turbulent channel, pipe, and boundary-layer 
flows for a wide range of Reynolds numbers. ^ 

We recall the boundary-layer approximation of the Leray-a model (see [7]). In the case of 
a zero-pressure gradient, consider a two-dimensional flow across a flat surface. Let x be the 
coordinate along the surface, y the coordinate normal to the surface, and u = (u, v) the velocity of 
the flow. Assuming that a is a function of x only, normalizing variables, and neglecting terms that 
are small near the boundary (see [l6[]), we arrive at a Prandtl-like boundary-layer approximation of 
the 2D Leray-a model: 



tW 
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v{x,y) = - —u{x,z)dz 
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(2) 



w = u 
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where (m, f ) are the components of the averaged velocity, p the averaged pressure, and w = 

The physical (non-slip) boundary conditions are u{x, 0) = v{x, 0) = 0, and 
{u{x,y),v{x,y)) (1,0) as y ^ oo. Looking for self-similar solutions to this system of the 
form 

a{x)=pS{x), ^ 

with 5{x) = 

1 

m = h — h" . (3) 



y 



X, we reduce Q to the following generalization of the celebrated Blasius equation: 



m'" + —hm" = 0, 



The physical boundary conditions for ([3]) are h{0) = h'(0) = and h'(^) — > as ^ — oo. 
This equation describes horizontal velocity profiles {h'(-)} in transitional and turbulent boundary 
layers with zero pressure gradients. In [[g] it was proved that the above boundary value problem 
has a two parameter family of solutions with the parameters being 

a := h"iO), b := h"'{0). 

From the derivation of ([3]) it follows that the averaged horizontal velocity profiles u for a fixed 
horizontal coordinate xo is modeled by 

u{xo,y) = ^h'(-y=], (4) 



\PVU 

where y is the vertical coordinate, Ue is the horizontal velocity of the external flow, his a solution 
to (|3]), P = (lim^^oo h'{y)Y^^, I is a local length scale that has to be determined, Ri = uj/u, and 
le is the external length scale le = v/ue- 

We now normalize quantities into wall coordinates 

y = — 5 u = — 

where 



/i - / ^ 
VP \l dy y=o 

and r is the shear stress at the wall. Writing dH) in wall units, a three-parameter family of velocity 
profiles is obtained ^ ;(■)• 



<uiy^) = ^h'(y^]. (5) 



III. SKIN-FRICTION CORRELATION 



A critical step in deriving a skin-friction correlation is reducing a three-parameter family of the 
velocity profiles ([5]) to a two-parameter family of turbulent velocity profiles, which is achieved 
with a help of the von Karman log law. 

For smooth surfaces, the mean velocity profile for the inner region is commonly approximated 
with the von Karman log law: 

= —Iny^ + B, 

K 
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where the von Karman constant, k ^ 0.4, and B ^ 5, are empirically determined constants. 
In the rough-wall case, the effects of uniform roughness are confined to the inner region, and 
are accounted for by modifying the semi-logarithmic part of the mean velocity profile. More 
specifically, Clauser [Si |91] showed that the semi-logarithmic region is displaced downward by 
an amount Au/ur- This amount of downward shift is commonly referred to as the roughness 
function, and represents the velocity defect from the standard velocity distribution over a smooth 
wall, and indicates the additional wall shear stress due to the roughness. Accounting for the 
roughness effect, the log law can then be written as 



u"*" = — In 

K 



{y + £)u^ 



B - AB, (6) 

V 

where e is the shift at the origin for the rough wall, y is measured from the top of the roughness 
element, and AB = Au/ur- The values of e and AB are determined by matching experimental 
velocity profiles with Q. 

We obtain turbulent velocity profiles by subjecting profiles © to three conditions of a weaker 
formulation of the von Karman log law: 

1. A turbulent velocity profile u^{y'^) has 3 inflection points in logarithmic coordinates. 

2. The middle inflection point oiu^{y^) lies on the line 

u+ = -\ny+ + B - AB. (7) 

K 

3. The line (|7]) is tangent to u^i^y^) at the middle inflection point. 
Let us fix a Reynolds number based on momentum thickness 

Re= r u+(l-^^\ dy+. (8) 



JO \ M^IOOJ/ 

Then dS]) and conditions (ii) and (iii) determine all three parameters a, 6, and / in dS]). Therefore, 
conditions (ii) and (iii) reduce ([5]) down to a two-parameter family of turbulent profiles {m^^ ^^}. 
Fig. [U shows the velocity profile for Rq = 700, AB = 1.57 and corresponding experimental data 
of Osaka et al. Jisl for d-type roughness. 

Note that the skin-friction coefficient is uniquely determined by a velocity profile cj = 
ab(^)^- Therefore, the skin-friction coefficient is now a function of Re and AB: 

Cf = f{Re,AB). (9) 
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FIG. 1: Velocity profile for Re = 700, AS = 1.57. 
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FIG. 2: Theoretical skin-friction correlation. 



rV. THE ROUGHNESS FUNCTION AND TRANSITION TO TURBULENCE 



The skin-friction law ^ is shown in Fig. [2] for several different values of roughness function, 
Ai? = Au/ur- At the critical point where Rg is at the minimum and c/ is at the maximal value, 
the second and third inflection points of the velocity profile collide, and the profile will then only 
have one inflection point for Rg < i?^"*. Therefore, the model indicates the minimal value of Rg, 
i.e. Rg^^^, for which a velocity profile can still be turbulent, i.e. the condition (i) is satisfied. One 
can interpret Rg^^^ as the minimal Reynolds number where the transition to turbulence may occur, 
often called the transitional Reynolds number in the literature. 

Fig. |3] shows how this critical Reynolds number depends on the roughness function, which was 
obtained by computing i?^"* for sixteen different values of AB. The following linear correlation 



6 



500 
400 
300 
200 
100 




I 1 1 1 1 1 r 



-51 .848AB+365.52 



J I I I I I I L 



0.1 



0.01 



-3 -2 -1 1 2 3 4 5 6 
AB 



0.001 



T 1 1 1 1 1 1 r 

0.00636°-^^^^ 



+ 

J----' 



J I I I I I I L_ 



-3 -2 -1 1 2 3 4 5 6 
AB 



FIG. 3: Critical value of the Reynolds number. FIG. 4: Maximum of the skin-friction coefficient, 
was found: 

Rf^ = -51.8AS + 365.5. (10) 

Let us now denote cf^^ = /(-R"'\ AB). Note that cj^^ is the largest possible value of the 
skin-friction coefficient in the turbulent boundary layer. As demonstrated in Fig. [2l higher val- 
ues of roughness function will allow velocity profiles to remain turbulent at lower Reynolds 
numbers, but furthermore will result in higher values of cj^"^. Fig. HI obtained by evaluating 
^max _ fi^Rcr\t^ AB), shows how the roughness function is predicted to influence the maximal 
value of the skin-friction coefficient. The following correlation was found: 



V. CONCLUSION 



Based on the Leray-a model of fluid turbulence, a generalized Blasius equation was formulated 
to describe streamwise velocity profiles in turbulent boundary layers with zero pressure gradients. 
Solutions of this fifth-order differential equation satisfying a weak formulation of the von Karman 
log law form a two-parameter family. The two parameters are the Reynolds number based on 
momentum thickness Re and the roughness function AB. This leads to a skin friction correlation 
Cf{Re, AB) and predictions of the critical Reynolds number Rg^"^{AB) and the maximal value of 
the skin-friction coefficient cJ^'^^AB). The critical Reynolds number is the minimal value of the 
transitional Reynolds number (which also depends on the intensity of the free-stream turbulence). 
In particular, it was shown that the greater AB is, the earlier the transition to turbulence may occur, 
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and furthermore the higher the skin-friction coefficient may peak. 
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